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In this Letter, a modified Chaplygin gas (MCG) model of unifying dark energy and dark matter 
with the exotic equation of state vmcg = BpMCG — cA — is constrained from recently observed data: 

T-H ' p MCG 

the 182 Gold SNe la, the 3-year WMAP and the SDSS baryon acoustic peak. It is shown that the 
best fit value of the three parameters (B,B s ,a) in MCG model are (-0.085,0.822,1.724). Furthermore, 
O |' we find the best fit w(z) crosses -1 in the past and the present best fit value w(0) = —1.114 < — 1, 

and the la confidence level of w(0) is -0.946 < iu(0) < -1.282. Finally, we find that the MCG 
model has the smallest Xmin value in all eight given models. According to the Alaike Information 
Criterion (AIC) of model selection, we conclude that recent observational data support the MCG 
' ■ model as well as other popular models. 

u: PACS numbers: 98.80.-k 

i-C 1 Keywords: Modified Chaplygin gas (MCG); dark energy; Alaike Information Criterion (AIC). 
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I. INTRODUCTION 

The type la supernova (SNe la) explorations [if, the cosmic microwave background (CMB) results from WMAP 
Q observations, and surveys of galaxies Q all suggest that the universe is speeding up rather than slowing down. 
The accelerated expansion of the present universe is usually attributed to the fact that dark energy is an exotic 
component with negative pressure. Many kinds of dark energy models have already been constructed such as ACDM 
quintessence Q], phantom [(J, generalized Chaplygin gas (GCG) [7], quintom [8], holographic dark energy tj], 



, and 

so forth. 

On the other hand, to remove the dependence of special properties of extra energy components, a parameterized 

•i-H . 

equation of state (EOS) is assumed for dark energy. This is also commonly called the model-independent method. 
The parameterized EOS of dark energy which is popularly used in parameter best fit estimations, describes the 
possible evolution of dark energy. For example, w = uio=const [lo| . w(z) = wq + w\z [ill ]. w(z) = wq + ji^ 121. 



W(z) =W + J^§y2 
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13], w{z) = i±£ ^i+W+*) _ i (here X = A 1 (l + z) + A 2 (l + zf + (1 - O 0m -A t - A 2 )) 
The parameters wq, W\, or A±, A 2 are obtained by the best fit estimations from cosmic observational datasets. 

It is well known that the GCG model has been widely used to interpret the accelerating universe. In the GCG 
approach, dark energy and dark matter can be unified by using an exotic equation of state. Also, a Modified Chaplygin 
gas (MCG) as a extension of the generalized Chaplygin gas model has already been applied to describe the current 
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accelerating expansion of the universe 0- 

The constraint on parameter B in MCG model, i.e., the 
added parameter relative to GCG model, is discussed briefly by using the location of the peak of the CMB radiation 
spectrum in Ref. 19j. In this Letter, we study the constraints on the best fit parameters (B, B s ,a) and EOS in the 



MCG model from recently observed data: the latest observations of the 182 Gold type la Supernovae (SNe) 20(, the 
3-year WMAP CMB shift parameter 2l| and the baryon acoustic oscillation (BAO) peak from Sloan Digital Sky 
Surver (SDSS) [22j. The result of this study indicates that the best fit value of parameters (B, B Sl a) in MCG model 
are (-0.085,0.822,1.724). Furthermore, we find the best fit w(z) crosses -1 in the past and the present best fit value 
w(0) = —1.114 < —1, and the la confidence level of u>(0) is —0.946 < w(0) < —1.282. At last, because the emphasis 
of the ongoing and forthcoming research is shifting from estimating specific parameters of the cosmological model to 
model selection 0], it is interesting to estimate which model for an accelerating universe is distinguish by statistical 
analysis of observational datasets out of a large number of cosmological models. Therefore, by applying the recent 
observational data to the Alaike Information Criterion (AIC) of model selection, we compare the MCG model with 
other seven general cosmological models to see which model is better. It is found that the MCG model has almost the 
same support from the data as other popular models. In the Letter, we perform an estimation of model parameters 
using a standard minimization procedure based on the maximum likelihood method. 

The Letter is organized as follows. In section 2, the MCG model is introduced briefly. In section 3, the best fit 
value of parameters (B,B s ,a) in the MCG model are given from the recent observations of SNe la, CMB and BAO, 
and we present the evolution of the best fit of w(z) with la confidence level with respect to redshift z. The preferred 
cosmological model is discussed in section 4 according to the AIC. Section 5 is the conclusion. 



II. MODIFIED CHAPLYGIN GAS MODEL 



For the modified Chaplygin gas model, the energy density p and pressure p are related by the equation of state 





^4 

Pmcg — Bp M ca 3 , (1) 

Pmcg 

where A, B, and a are parameters in the model. 

Considering the FRW cosmology, by using the energy conservation equation: d(pa 3 ) = — pd(a 3 ), the energy density 
of MCG can be derived as [18j 

Pmcg = Pomcg[B s + (1 - B.)(l + z fV+ B ^+«)]Th, (2) 
for A ^ -1, where a is the scale factor, B s = A i+a . In order to unify dark matter and dark energy for the 

{1+B)p 

MCG model, the MCG fluid is decomposed into two components: the dark energy component and the dark matter 
component, i.e., pmcg = Pie + Pdmi Pmcg — Pde- Then according to the relation between the density of dark matter 
and redshift: 



Pdm — P0dm(l + z) , 



(3) 



the energy density of the dark energy in the MCG model can be given by 

Pde = Pmcg - Pdm = Pomcg[B s + (1 - + zf^+^+^rh - Podm (l + z f. (4) 

Next, we assume the universe is filled with two components, one is the MCG component, and the other is baryon 
matter component, ie., pt = Pmcg + Pb- The equation of state of dark energy can be derived as 

_ (1 - n Qb )[B s + (1 - B s )(l + z )3(i+B)(i+a)]-Tfe[- jB3 + B (l - B s )(l + z )3(i+B)(i+q)] 
yj(i e — 1 , (o) 

(1 - n ob )[B s + (1 - + z )3(l+B)(l + a)]— _ nodm{1 + z) 3 

where £lodm and f^ofc are present values of the dimensionless dark matter density and baryon matter component. 
Furthermore, in a flat universe, making use of the Friedmann equation, the Hubble parameter H can be written as 

H 2 = ^=H 2 E\ (6) 

where E 2 = (l-Q. oh )[B a + (1 - B s )(l + ^) 3 ( 1 + B )( 1 + Q; )] tts + Q, ob (l + zf. H denotes the present value of the Hubble 
parameter. When B = 0, equation ^ is reduced to the GCG scenario. 

In the following section, on the basis of equation ([S]), we will apply the recently observed data to find the best fit 
parameters (fiobi B, B s ,a) in MCG model. For simplicity, we will displace parameters (Clob, B, B s ,a) with 6 in the 
following section. 



III. THE BEST FIT PARAMETERS FROM PRESENT COSMOLOGICAL OBSERVATIONS 

Since type la Supernovae behave as Excellent Standard Candles, they can be used to directly measure the 
expansion rate of the universe up to high redshifts (z > 1) for comparison with the present rate. Therefore, they 
provide direct information on the universes acceleration and constrain the dark energy model. Theoretical dark 
energy model parameters are determined by minimizing the quantity 

2 itt a\ (Pobs(zi) - pth(zi)) 2 , . 

XSNei H 0,O) = ~2 > ( 7 ) 



where N = 182 for the Gold SNe la data 0, 

a 1bs-t are errors due to flux uncertainties, intrinsic dispersion of SNe la 

absolute magnitude and peculiar velocity dispersion respectively. The theoretical distance modulus p t h is defined as 

Hth(Zi) = mth(zi) — M = hlog w {D L (z)) + Uog w {-^—) + 25, (8) 

Mpc 

where 

D L (z) = H 0dL (z) = (l + Z )l ^^-y (9) 

Pobs is given by supernovae dataset, and is the luminosity distance. 

The structure of the anisotropics of the cosmic microwave background radiation depends on two eras in cosmology, 
i.e., last scattering and today. They can also be applied to limit the model parameters of dark energy by using the 



i—\ ® ~obs;z 



4 



shift parameter 24 [ 



R = v 7 ^ 



Om 



Hr\d,Z 



(10) 



;o H(z';H ,9Y 

where z rec = 1089 is the redshift of recombination, fiom is present value of the dimensionless matter density, inclu ding 
dark matter and the baryon matter component. By using the three-year WMAP data 



R = 1.71 ±0.03. 



25] , R can be obtained as [26 1 

(11) 



From the CMB constraint, the best fit value of parameters in the dark energy models can be determined by minimizing 

(R(H ,e) - 1.71) 2 



0.03 2 



(12) 



Because the universe has a fraction of baryons, the acoustic oscillations in the relativistic plasma would be imprinted 
onto the late-time power spectrum of the nonrelativistic matter 27J. Therefore, the acoustic signatures in the large- 



scale clustering of galaxies can also serve as a test to constrain models of dark energy with detection of a peak in the 
correlation function of luminous red galaxies in the SDSS . By using the equation 

HqcIz 



A = \/£la m E{zBAo) 



-l/3f 



1 



12/3 



zbao Jo H(z';H ,6) 

and A = 0.469 ± 0.017 measured from the SDSS data, zbao = 0.35, we can minimize the x^Qdefmed as [2 

2 (T r n (A(z';H Q ,0)-OA69) 2 
Xbao{Ho,V) ixoifi • 

As one can find that the gravitational clustering in MCG model presented in Rcf. 



(13) 



(14) 



Uj , ensures that the observational 



datasets from CMB and BAO can be applied to constrain the MCG model. Hence, we combine these three datasets 
to minimize the total likelihood Xtotai 



Xtotal 

(H , 6) - XsNe + XCMB + XbAO- 



(15) 



On the one hand, since we are interested in the model parameters 0, the Hq contained in Xtotai (Ho, 0) is a nuisance 
parameter and will be marginalized by integrating the likelihood L(6) — J dHoP(Ho) exp (— x 2 (Ho, 0)/2), where 
P(Hp) is the prior distribution function of the present Hubble constant, and a Gaussian prior Ho — 72±8kmS~ 1 Mpc -1 
[29j is adopted in the Letter. We know that the prior knowledge of cosmological parameter fiob has been obtained 
by several other observations, such as fiofc/i 2 = 0.0214 ± 0.0020 from the observation of the deuterium to hydrogen 
ratio towards QSO absorption systems 



flobh 2 



0.021 ± 0.003 from the BOOMERANG data 



31) and Cl 0b h 2 



32j for the observation of CMB. Thus, in order to get the interesting result for the 



0.022±{J$g from the DASI results 

value of f2ob) we treat flp h as a free parameter with Gaussian prior distribution centered in rip" 6 with spread ctq 
And following Ref. 33J, the "weak" prior for par ameter fiob will be used in our analysis, i.e., let it have a relative 
larg er variable range ttobh 2 



0.0214 ±0.0060 



Xtotal-prior(@) ~ Xtotal(@) ± 



Thus, the Xtotai(Ho,6) m equation (|T5|) will be reconstructed as 
(ftob - Slpb" 6 ) 2 



(T, 



(16) 



Qob —prior 



where Xtotai(@) denotes the total x 2 (#) obtained without imposing prior knowledge of fl 



06 • 
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FIG. 1: The best fits of w(z) with la confidence level (shaded region) 



By using the maximum likelihood method for equation (|16p. we obtain the best fit values (f^ofcj B, B Sl a) in the 
MCG model (0.041,-0.085,0.822,1.724) with X 2 min = 157.272. Figure 1 shows the la confidence level of the best 
fit w{z) calculated by using the covariance matrix. From Fig.l, it is easy to see that the best fit w(z) cross -1 at 
about z — 0.140 and the present best fit value w(0) = —1.114 < —1. Obviously, it can be shown that the fact 
that w(z) cross over the boundary of w = — 1 in MCG model is consistent with the results given by Refs. [l^l [3 
j^ , where w(z) crossing -1 is first found using SNe data. Furthermore, we obtain the la confidence level of w(0), 
—0.946 < w(0) < —1.282. The possibility of w(0) > —1 camrt be excluded in la level. At last, it can be seen that 
the cosmological constant model (i.e.,w(z) = —1) is not in la confidence contour of the best fit dynamical w(z). 



IV. THE PREFERRED COSMOLOGICAL MODEL 



It is interesting to ask which model of an accelerating universe is preferred by recently observed data over 
many models. We also want to know how well the MCG model fits the recently observed datasets as compared to 
other models. We make use of the values of Xmin an d the objective Alaike Information Criterion (AIC) to solve the 
questions above. 



On the basis of the description in section 3, we obtain the values of Xmin by minimizing the Xt 



otal— prior 



the 



corresponding models, where ftom is treated as a free parameter with a Gaussian prior in the range fiom — 0.29 ±0.07 
Q for the ACDM model and model-independent cases, and the results are listed in Table 1. It can be seen that the 
MCG model has the smallest Xmin value. Table 1 contains the best fit parameters corresponding to the different 
models. 

In cosmology the AIC was first used by Liddle [3^], and then in subsequent papers [38 



. It is defined as 



where C r 



AIC = -2 In C{6 I data) max + 2K, (17) 
is the highest likelihood in the model with the best fit parameters 9, K is the number of estimable 
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case model 


/Cmin 


Best fit parameters 


ACDM 


162.302 


fl 0m = 0.286 


w=constant 


159.962 


Q 0m = 0.288, w = -0.870 


w(z) — WQ + W\Z 


159.765 


Ho™ = 0.292, w = -0.893, wi = 0.009 


w(z) = w + ^ 


158.635 


tl 0m = 0.289, w = -1.041, wi = 0.751 


w(z) - w + {1+ \ y2 


157.715 


Ho™ = 0.282, w = -1.314, wi = 3.059 


w{z) - 1+ 3 * dl±^kii±£) 1 


159.067 


n 0m = 0.292, Ai = -0.302, A 2 = 0.188 


GCG 


159.444 


f2 6 = 0.041, As = 0.678, a = -0.136 


MCG 


157.276 


Slob = 0.041, B = -0.085, B s = 0.822, a = 1.724 



TABLE I: The values of Xmini an d best fit model parameters against the model 



parameters (0) in the model. The term — 21n£(0 | data) in Eq.(fl~7f is called % 2 and it measures the quality of model 
fit, while the term 2K in Eq. ([T7l) interprets model complexity. For more details about AIC please see Refs. H 

m H • 

In what follows, we will estimate which model is the better one for all the models in Table 2. The value of AIC has 
no meaning by itself for a single model and only the relative value between different models are physically interesting. 
Therefore, by comparing several models the one which minimizes the AIC is usually considered the best, and denoted 
by AIC m i n =min{ AIQ, i — 1,...,N}, where i — 1,...,N is a set of alternative candidate models. The relative 
strength of evidence for each model can be obtained by calculating the likelihood of the model £(M, | data) oc exp 
(— Aj /2), where Aj = AIQ— AIC m i„ over the whole range of alternative models. The Akaike weight Wi is calculated 
by normalizing the relative likelihood to unity and corresponds to posterior probability of a model. The evidence for 
the models can also be judged by the relative evidence ratio ^ = £(Mj | data) • ^ model i is the best one, the relative 
evidence ratio gives the odds against the model. The rules for judging the AIC model selections are as follows: when 
< Ai < 2 model i has almost the same support from the data as the best model, for 2 < Aj < 4, model i is 
supported considerably less and with Aj > 10 model i is practically irrelevant. 



case model 


AIC 


Aj 


Wi 


Odds 


ACDM 


164.302 


0.587 


0.149 


1.342 


w— const ant 


163.962 


0.247 


0.177 


1.130 


w(z) = Wo + WlZ 


165.765 


2.050 


0.072 


2.778 


l \ | Till Z 

w(z) — w + 


164.635 


0.920 


0.126 


1.587 


W{Z) =W + JYf^yJ 


163.715 





0.200 


1 


w(z) - *+* ^+' A / 1 +') l 


165.067 


1.352 


0.102 


1.961 


GCG 


165.444 


1.729 


0.084 


2.381 


MCG 


165.276 


1.561 


0.091 


2.198 



TABLE II: The value of AIC, Akaike difference, Akaike weights Wj and odds against the model 



Thus based on the values of Xmin °f a ^ models in Table 1, the evidence of the AIC can be calculated. We find that 
the best model is the one following w(z) = wo + r^t§^ i n terms of its AIC value. Taking it as a reference, we calculate 
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the differences between the models by using the AIC differences Aj, Akaike weights Wi and odds against alternative 
models. Table 2 gives the calculating results. Note that the model selection provides quantitative information to 
judge the "strength of evidence", not just a way to select only one model. From Table 2 it is easy to see that, the 
MCG model has almost the same support from the data as the best model, because the value of A$ for it is in the 
range 0-2. Furthermore, it can be shown that the recent observational data supports all of the models in Table 2 
except for the case of w(z) — wq + wiz since the value of Aj for this case is a little bigger than 2. It has a less support 
from recent observations. On the other hand, we can see that the MCG model is favored by observational data more 
than GCG model according to the Akaike weights Wi in Table 2. Finally, the odds indicates the difference between 
MCG model and the best one is 2.198 to 1. 



V. CONCLUSION 



In summary, the constraints on the MCG model, proposed as a candidate of the unified dark matter-dark energy 
scenario, has been studied in this Letter. We obtained the best fit value of the three parameters (B,B s ,a) in the 
MCG model (-0.085,0.822,1.724). Meanwhile, it is easy to see that the best fit w(z) can cross -1 as it evolves with 
the redshift z, and the present best fit value w(Q) = —1.114 < —1. Furthermore, it is shown that the la confidence 
level of w(0) is —0.946 < w(Q) < —1.282, and the possibility of w(0) > —1 camrt be excluded in la level. We can 
see that the cosmological constant model (i.e.,w(z) — —1) is not in la confidence contour of the best fit dynamical 
w(z). Finally, in order to find the status of MCG scenario in a large number of cosmological models, we compared 
the MCG model with other seven popular ones offering explanation of current acceleration of the universe in terms of 
the values of Xmin an d AIC quantity. We find that, as the quantity Xmin measures the quality of model fit, the MCG 
model is preferred by recent observational data because of its a small minimum \ 2 value. On the other hand, it is 
shown that the MCG model has a slightly high value of AIC due to its many parameters. However, according to the 
rules of judgment of the AIC model selection, we conclude that recently observed data supports the MCG model as 
well as other popular models, because the value of Ai for it is in the range 0-2 relative to the best model. In addition, 
the result of study shows that the recent observational data equivalently supports all of the models in Table 2 except 
for the case of w(z) = wq + w\z. We expect the new probers such as SNAP and Planck surveyor can provide more 
accurate data and further explore the nature of dark energy. 
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